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It  is  well  knownthat  moment  estimates  for  the  Weibull  distribution  with 
shape  parameters  $  and  scale  parameter  6  can  be  found  by  solving  for  B  the 
equation 

s2/x2  =  Er(l  ♦  2/B)/r2(l  +  1/g)]  -  l;  (1) 
and  having  found  6,  obtaining  the  estimate  for  6  as 
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<S  -  x/r(l  +  1/8),  (2) 

where  x  and  s  are  the  sample  mean  and  variance. 

To  facilitate  the  solutions  for  6  and  6,  we  provide  a  table  of  g  vs. 
the  right-hand  side  of  equation  (1)  and,  in  parallel,  of  r(l  +  1/B). 

The  method -of -moments  is  loiown  to  be  inefficient  with  respect  to  a 
number  of  other  estimation  procedures  and  should  be  used  only  when  the 
nature  of  the  data  available  does  not  permit  the  use  of  one  of  these 
better  methods,  or  if  the  sample  size  is  "large".  Examples  of  the  former 
situation,  and  of  the  use  of  the  tables,  are  provided. 
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W.R.  Blischke,  L.  Guin,“  E.M.  Scheuer’ 

1.  Introduction  and  Summary 

It  is  well-known  that  moment  estimates  for  the  Weihull  distribution 
with  shape  parameter  B  and  scale  parameter  <S , 

B-l 

(IJ  f(t)  =  j-  (j)  exP  ,  B  >  0,  6  >  0 

can  be  found  by  solving  for  6  the  equation 

(21  82/02  -  [  r (i  +  2/g) / r2(i  +  i/B)  ]  -  i 

and,  having  found  0,  obtaining  the  estimate  for  6  as 

(3)  ^  =  y/rfl  >  1/BO • 

Inc  quantities  y  and  a2  are  sample  estimates  of  population  mean  and 
variance. 

To  facilitate  the  solution  for  B  and  6  we  provide  herein  a  table  of 
f->  vs.  the  right-hand  side  of  equation  (2)  and,  in  parallel,  of  T(1  +  1/B). 

2 .  Why  use  the  method-of-moments? 

It  is  well-known  that  the  method-of-moments  is  inefficient  with  respect 
to  a  number  of  other  estimation  procedures  and  that  it  should  be  used  only 
when  the  nature  of  the  data  available  docs  not  permit  the  use  of  one  of 
these  better  methods,  or  if  the  sample  size  is  "large."  Here  is  an  in¬ 
stance  of  the  former  situation. 

The  data  available  for  analysis  come  from  k  independent  samples,  all 
presumably  from  the  same  Weihull  population.  The  sample  sizes  are 
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m ,  n2 ,  n^.  The  individual  values  Ty  (j  =  1,2,  ...,  n2 ;  i=l,  k) 

ni 

are  not  available;  only  the  sums  ^  Tjj  --  representing  the  total 

j-1 

time-on-test  in  the  i-th  sample,  and  the  sample  sizes,  nj  (i=l,2 . k) 

are  known. 

k  k 

Let  T  =  ^  Tj  and  n  =  ^  n-j.  It  is  easy  to  verify  that  yj  --  Tj/rij 
1  1 

is  an  unbiased  estimate  of  y,  the  population  mean,  based  on  data  from  the 
i-th  sample  alone,  that  y  =  T/n  is  an  unbiased  estimate  of  y,  and  that 

k 

(4)  a2  =  £7y  ni(yi  '  k1)2 
1 

is  an  unbiased  estimate  of  o?,  the  population  variance.* 

An  example  follows: 

.i.  example 

We  generated  random  variates  from  a  Weibull  distribution  with  6  =  1000 
and  8=3.  We  took  k  =  4,  n*  =  5,  n2  =  8,  n5  =  9,  n»  =  5.  We  found 
Tt  =  5525.22,  T2  =  5321.20,  T3  =  10,783.55  and  T„  =  5806.61.  Thus  y, 

=  1105.04,  y?  =  665.16,  y3  =  1198.17,  y„  =  1161.32,  y  =  1016.17, 

.'2  =  4’26,213.29  and  o2/y2  -  .461177.  Using  our  table,  wc  obtain 

,  t  1016.17  _  ,nr  ,  , 

;■  ~  1 .  :a,  <uiJ  <5  =  ~~ 90275"  ~  1125.  ('4. 

4 .  The  Tables 

As  stated  at  the  outset,  we  tabulate  8  vs  £r(l  +  2/B j/T2  f 1  +  1/8 ij  -1 
vs  r(l  +  1/P)  for  P  over  the  range  0.200  to  10.000,  this  chosen  because 

*lr  should  lie  noted  that  the  foregoing  doesn’t  depend  at  all  upon  the  Wcihull 
assumption.  It  is  true  generally  under  the  assumptions  of  independence  and 
homogeneity  -  provided,  of  course,  that  the  common  distribution  has  a  variance. 


it  covers  what  we  believe  to  be  the  set  of  values  likely  to  be  encountered 
in  practice.  The  increments  for  0  are  not  uniform,  but  selected  to  provide 
for  "smooth"  coverage  of  the  functions  being  tabulated.  The  tabulation  is 
l'or  $  =  0.2(. 001)0. S(. 01)1. 0(. 10)10.0. 

The  tables  are  appended. 

S.  Other  Treatments  of  This  Problem 

Marquina  (1979)  has  suggested  a  root-finding  procedure  to  estimate  e 
from  eq.  (2)  --  or  rather  its  reciprocal.  Shooman  (1968)  has  suggested 
an  iterative  approach.  The  latter  seems  less  tractable  than  our  tabular 
approach.  If  more  accuracy  is  needed,  our  approach  could  be  used  to 
obtain  a  good  starting  value  for  a  root-finding  procedure,  such  as  the 
one  proposed  by  Marquina. 
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